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Abstract 

We investigate the production of Higgs bosons A/h via the Yukawa process 
Z — > ff — > ffA/h taking into account QCD corrections and fermion mass effects. 
We estimate the discovery reach of this process and/or bounds on the parame- 
ter space of two-Higgs doublet models which can be derived from the analysis of 
bbr + T~ events at LEP1. These limits have important consequences for extended 
models of electroweak interactions and their experimental verification/improvement 
are welcome. 



1 Introduction 



Experimental searches at LEP1 for the neutral Higgs bosons in the context of extended 
models of electroweak interactions have been performed in the Bjorken [Z — > Z*h) and 
pair production (Z — > Ah) processes. However, the Yukawa process (Z — > ff — > ffA/h), 
in which the Higgs boson is radiated off the final fermion, provides an alternative Higgs 
production mechanism. Although the possible importance of this process in the context 
of two-Higgs doublet models, both in supersymmetric and non-supersymmetric versions, 
has been pointed out long time ago P], 0, [3| , experimentally it has not been looked for so 
far. 

Recently there has been a renewed interest in the Yukawa process in the context of 
attempts to explain experimental measurements of Rb, the branching ratio for Z decays 
to bb II, ||. In particular, in the paper [Q it has been argued, based on (among others) a 
bound on parameter space tan/3 — Ma derived from the Yukawa process Z — > bbA followed 

1 Supported in part by Polish Committee for Scientific Research and the EC grant under the contract 
CHRX-CT92-0004 
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by A — ► bb [|| , that problem of Rb in the supersymmetric model could be explained only 
for tan/3 ~ 1. This bound has been obtained from the simple assumption that if 20-30 
bbA events had been produced at LEP1 it is likely to have been noticed because the 
QCD background from Z — > 46 is small. However, the analyses of 46-jet events, for 
example in Z — > hA — > 4b searches §, show that the background from Z — > bbgg due to 
misidentification of 6-quark jets is much more important. Although the Yukawa process 
has a different topology, similar background problems may appear in its experimental 
analysis and a definite conclusion can be drawn after a dedicated experimental analysis 
in 46 channel will have been performed. On the other hand, the bbr + r~ channel is almost 
background free and therefore it can be exploited to estimate limits on the parameters of 
the extension of the Standard Model. 

In this paper we estimate the parameter bounds and/or discovery reach of the Yukawa 
process using the bbr + r" channel. To this end we provide the simple formulae including 
all fermion mass terms for the matrix element for process e + e~ — > f ftp for both (f) = h 
and (f) = A. They may prove useful in constructing Monte Carlo programs in order to 
include experimental constraints, efficiencies etc. 

In the numerical analysis we also include QCD corrections which turn out to be very 
important. Assuming that no signal events are found we obtain weaker bounds on tan/3 — 

than in || . Nevertheless our results demonstrate the physics potential of the Yukawa 
process: even if the Higgs boson is not found in this process, the exclusion plot derived 
from negative search will have interesting consequences for extensions of the Standard 
Model. 



2 Higgs particle production at LEP1 

In the minimal extensions of the Standard Model (SM), either supersymmetric (MSSM) 
or non-supersymmetric (2HDM), two Higgs doublets are employed 0]. In this paper we 
consider 2HDM of type II in which, like in the MSSM, one of the Higgs doublets couples 
to up-type, and the other to down-type fermions. After the symmetry breakdown there 
are five physical Higgs particles: two CP-even neutral scalar Higgs bosons denoted by h 
and H (we assume that M h < M H ), a CP-odd neutral - A (usually called a pseudoscalar) , 
and a pair of charged particles H ± . Their experimental discovery and measurement of 
their couplings provides one of the major motivations for and constitues a primary physics 
goal to be achieved at present and future colliders. 

Currently the best experimental limits on Higgs boson masses are obtained from LEP1 
experiments performed at the Z boson peak. At the tree level the neutral Higgs bosons 
h and A can be produced in Z boson decays via three processes: 

Z Z*h (1) 
Z ^ Ah (2) 
Z - ff^ffA/h (3) 

In the SM the associated Higgs pair production process (fj) is absent and the Yukawa pro- 
cess (||) is suppressed by a factor (rrif/M w ) 2 with respect to Higgs-strahlung process(|I]). 
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Therefore the lower limit on the Higgs mass Mh sm > 65.2 GeV has been derived || from 
the negative search via the process (0) at LEP1. 

In the two-doublet extensions of SM the Higgs-strahlung process ([!]) occurs at a lower 
rate 

T(Z^hZ*) = T S M(Z^H SM Z*)sm 2 ((3-a) (4) 

because the ZZh coupling is reduced by a factor sin 2 (/3 — a), where a and (3 are the 
mixing angles in the neutral and charged Higgs sectors, respectively. The negative search 
in this process results in the upper limit |J on sin 2 (/3 — a) for a given Mh, for example 
sin 2 (/? — a) < 0.1 if M/, < 50 GeV. However, for small sin 2 (/? — a) one can exploit the 
Higgs pair production process (0), if kinematically allowed, because it is complementary 
to the process ([I]) in the sense that 

T(Z -> hA) = 0.5T(Z -> vu) cos 2 (/3 - a)A 3/2 , (5) 

where A = (1 — Kh — ka) 2 — ^k\k 2 a , with «j = mf/Mf. Such a procedure has been adopted 
at LEP1. An independent contraint on the non-SM contributions to Z decay is obtained 
from the analysis of the Z-lineshape, where the process (5) also contributes. 

In this paper we consider a scenario in which the process (|1|) is suppressed (small 
sin 2 (/3 — a)) and the process (0) is closed kinematically {m h + rriA > Trig)- Then the 
Yukawa process (3) becomes a dominant Higgs production mechanism at LEP1. Such a 
scenario can easily be realized in the 2HDM where the Higgs masses and mixing angles are 
not constrained theoretically. On the other hand, in the MSSM, due to supersymmetry 
relations, the Yukawa process may become dominant only for > 50 GeV and tan (3 > 10 

!■ 

The number of expected events crucially depends on hff and Aff couplings: the SM 
Higgs coupling {\/2GF) l ^ 2 mf is multiplied by the factors <fy// given in Table 1, where D 
(U) is a generic notation for down-type (up-type) fermions, respectively. 



9<t,ff 


f = D 


f = u 


<p = h 
(f> = A 


— sin a/ cos (3 
tan p 


cos a/ sin P 
l/tan/3 



Table 1: Higgs boson couplings to fermions relative to the SM Higgs couplings. 

In the interesting large tan P case the couplings of both A and h to down-type fermions 
are strongly enhanced by the same factor tan/3 because small sin 2 (/? — a) implies^, to a 
good approximation, that sin a/ cos/? ~ tan/?. In such a case the Yukawa process can be 
a source of a significant number of Higgs bosons h or A, whichever is lighter. For example, 
at LEP 1 one can expect fllj for M^/a = 10 GeV and tan/3 = 20 a few thousand of bbh or 
bbA events in 10 7 Z decays which corresponds to the Z partial decay width of ~ 1 MeV. 

2 A detailed discussion will be given elsewhere . In what follows we take a = (3 for the hbb and Iitt 
couplings. 
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3 The Yukawa process 



The rate for the process is very sensitive to the fermion mass effects both in matrix 



elements and in the kinematical limits as demonstrated in our previous paper fLl| . The 
analytical formulae for the matrix element of the process e + e~ — > ff<p, including all 
fermion mass terms, can be written in a compact form. We present them in Appendix, as 
they may appear to be useful in constructing Monte Carlo programms in order to include 
experimental effects. 

In this paper we consider the effect of QCD corrections in the Yukawa process. In the 
case when the Higgs boson is radiated off the b quark (/ = b) and in hadronic decays 
(0 — > bb) the QCD corrections are included by using the running quark mass rrif = m^M^) 
in the Higgs-quark coupling. This has an important effect on the production cross section. 
For example, for = 50 GeV the cross section is reduced by a factor 2 in comparison 
to fixed quark mass case, i.e. neglecting QCD corrections. 

We observe that the cross section for the Yukawa process is proportional to tan 2 (3. 
The additional dependence on tan/3 that appears in the branching ratios for dominant 
Higgs decay modes (f) — > t + t~ and bb for tan/3 > 5 is very weak and therefore can be 
neglected^. To derive the limits on tan (3 — (with <p = h,A) we proceed as follows: 

a) We perform the calculations at the Z pole for 10 7 Z's produced. 

b) We consider two different final state configurations: bb(r + r~) and T + r~(bb), where 
the first pair denotes the fermions produced in e + e~ collisions which radiate the 
Higgs boson (either h or A) and the pair in the parenthesis is coming from Higgs 
decay. 

c) As we discussed, the bbr + T~ channel seems to be the most promising both from the- 
oretical (negligible QCD background) and experimental (relatively good efficiency) 
points of view. For example, the L3 Collaboration [I3[] quotes the acceptance for 



signal events from Z — > hA — > bbr + r of the order of 5%. We assume that in the 
analysis of the Yukawa process similar efficiency can be achieved 

d) We assume no signal events and we derive the bounds corresponding to 95% CL. 

The results are presented in Fig.l for the pseudoscalar Higgs boson A and in Fig. 2 
for the scalar h. Parameters above the lines are excluded. The solid lines represent the 
limits obtained with running 6-quark mass mb(M<p), and the dashed lines correspond to 
the case of fixed 6-quark mass m& = 4.6 GeV, i.e. neglecting QCD corrections. The 
QCD corrections lower the Higgs-quark coupling and consequently weaker bounds can be 
established. For comparison the dotted lines delineate the excluded region based on 20 
46-jet events expected. 

Note that in the MSSM case the mass range of interest for Mh and Ma is lying 
above, say 50 GeV. On the other hand in the 2HDM relatively light Higgs bosons are 
still not excluded, for example for tan (3 = 20 the Higgs masses as low as 25 GeV, and for 
tan p = 55 masses of 45 GeV are not excluded. 



We use the program developed in ref.[12| to calculate the Higgs decay branching ratios. 
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To summarise, the Yukawa process provides interesting limits on the parameter space 
of 2HDM and MSSM. Both QCD and fermion mass effects influence the results consider- 
ably. 

Although the first attempt to include the Yukawa process in Higgs searches at LEP1 



has been made recently by L3 Collaboration ]X3|] nevertheless this process still awaits a 
serious experimental investigation. 
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Appendix 

We consider the Yukawa process 

with / = b or r and the corresponding 4-momenta are given in brackets. The electron 
mass is neglected (kf 2 = m>l = 0) and the final state fermion (pole) mass is denoted by 
m {p\ 2 = rn 2 ). The Higgs mass is denoted by and s — (ki + k 2 ) 2 - 
The differential cross section EiE 2 da/d 3 pi d s p 2 reads 

da 

''d 6 p^p 2 128tt 5 s 2 V 



^2 j3 _ j3 _. = mo 5 2^ \ M \ ~ 2 + Si + S 2 ) ( 6 ) 



_ 2-^1 _ 2-^2 _ 22^ 

X l — — 1=1 X 2 — — l=i x <h — 7=5 {<) 

VS y/S Y y/S 

where E^ 2 and denote CM energies of fermions and the Higgs boson 0, respectively. 

The matrix element, averaged over the initial and summed over final fermion polar- 
izations, can be cast in the following form 

2 G 
]M\ = Q^ 2 a 2 el N c -^gl ff -m 2 f s 

sn 2 n 2 2! + ( a e + vl)[vfr 2 + a /C 2 ] + 2a e v e a f v fX 

We^f ^ + ( s _ M 2 )2 + M 2p| 

, n n {Ivevn 2 + Q e a/X)(g ~ M 2 ) 

s{{s - M 2 ) 2 + M 2 T 2 Z ) h [ ) 

The color factor is N c = 3 for / = b and N c — 1 for f — r. We distinguish the fermion 
mass which enters the coupling, m/, from the pole mass m. In the case of the r lepton 
we take rrif = m = m T . However, in the case of the 6-quark, we include QCD corrections 
by taking the running quark mass mj = m^M^). Q e ,Qf are electric charges of electrons 
and fermions, respectively (in units of e) whereas a e (af),v e (v e ) correspond to axial and 
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vector couplings of electron (fermion /) to Z-boson with a = Izl /2 sin 6w cos 6w and 
v = (I 3L — 2Q sin 2 9 w )/2 sin 9 W cos 9 W . 

The contributions of the photon, Z-boson exchange and the 7 — Z interference are 
denoted by 7 2 , ( 2 and Xi respectively, and can be written as follows: 
a) Pseudo scalar 

In the case of the pseudoscalar production, = A, we have g\^j = tan 2 (3, = M A and 



7272 -%a — + ^( n i6 + n 2 6)] + t-t^a - 2 (— 1 



1 Mi 2 2m 



2 



s 



[4 ^(n^-n^i)] 



1 r 2 tt o/ M l- 4m 1 x 21 

+7^ [4i + + 2 ^ + 1 - x A )x 2 T ] 



1 ^2 ,2 - ,21 1 
t2t2~L"l?2 - -2?l] - 7^ 
■51^2 S U« 



X = 7272 — [S^-^l-TVK^-^O-^An] (9) 
b) Scalar 

For the scalar, = /i, we have g 2 ^ = (cos a/ sin/5) 2 , M<p = M h and 
o 1 M? — Am 2 r 9 —2m 2 — s . „ ^ 

7 = 7272^ [4 — +^(ni6 + n 2 6)] 

S1S2 6 6 

+7^[^-2(— + 1)4] 

C = TiTii— [4 ^(n^-n^i)] 

S1S2 S S 

[4^^ + ^(ni6 + n 2 6)]} 
+t^[4 + - 16— (x h — + 2(— ^ + 1 - x h )4] 

?1?2 S S S 

X = 7272 [— (Si^-^i 2 ) - S 2 Ci)] 

?1?2 s s 

' [(5x6 - s 2 eO - x ft An] (10) 



66 

where the following variables are introduced 



X21 = x 2 - x 1 (11) 

Sl = l-X!, 6 = 1-^2 (12) 

4 = 4(p! + p 2 ) ■ h ( Pl + p 2 ) ■ k 2 /s 2 -l + Xh - Mils (13) 

U 1 = 8(p 1 -k lPl -k 2 )/s 2 , U 2 = 8(p 2 -k lP2 -k 2 )/s 2 (14) 

n_ = rii-n 2 (15) 

Si = 4pi • (k± - k 2 )/s, E 2 = Ap 2 -(k 1 -k 2 )/s (16) 

All = 8(pi • k x p 2 ■ k 2 - pi ■ k 2 p 2 ■ ki)/s 2 (17) 
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Figure 1: The tan/3 — M A exclusion/discovery plot for the Yukawa process for A produc- 
tion. The region above the lines is excluded if no signal in bbr + r~ mode is observed. The 
solid line represents the results when QCD corrections are accounted for by taking the 
running 6-quark mass and the dashed line corresponds to fixed 6-quark mass m h = 4.6 
GeV. For comparison the dotted line shows the region that can be excluded based on 20 
events with 5-quark jets (using the running 6-quark mass). 
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Figure 2: The tan (3—M h exclusion/discovery plot for the Yukawa process for h production. 
The region above the lines is excluded if no signal in bbr + T~ mode is observed. The solid 
line represents the results when QCD corrections are accounted for by taking the running 
6-quark mass and the dashed line corresponds to fixed 5-quark mass mj = 4.6 GeV. For 
comparison the dotted line shows the region that can be excluded based on 20 events with 
6-quark jets (using the running 6-quark mass). 
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